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Introduction 

The  global  ocean  has  its  own  "weather"  phenome¬ 
na,  although  with  greatly  different  time  and  space 
scales  compared  to  the  atmosphere.  Oceanic  mesoscale 
eddies  are  typically  about  100  km  in  diameter  which 
makes  them  20-30  times  smaller  than  comparable 
atmospheric  highs  and  lows.  The  ocean's  "jet  streams" 
are  the  western  boundary  currents  and  their  extensions 
into  the  interior  ocean.  The  currents  have  speeds  on  the 
order  of  1  m/s  compared  fo  afmospheric  speeds  fhaf 
can  be  100  fimes  fhis  value.  The  space  scales  of  fhe 
meanders  on  fhese  high-speed  sfreams  are  similar  fo 
fhose  for  fhe  eddies  mentioned  above.  Knowing  and 
predicfing  fhese  oceanic  mesoscale  feafures  have 
numerous  naval  applicafions  which  include  facfical 
plarming,  optimum  frack  ship  routing,  search  and  res¬ 
cue  and  supplying  boundary  conditions  for  high  reso- 
lufion  coasfal  models,  fo  name  a  few. 

However,  operational  global  ocean  modeling  has 
lagged  far  behind  afmospheric  modeling  because  fhe 
smaller  space  scales  mean  approximafely  four  orders  of 
magnifude  more  compufer  fime  and  fhree  orders  of 
magnifude  more  compufer  memory  are  required.  In 
addition,  unlike  fhe  mefeorological  radiosonde  nef- 
work  fhaf  provides  inifial  condifions  from  fhe  surface 
fo  fhe  fop  of  fhe  afmosphere,  fhere  are  very  few  obser- 
vafions  below  fhe  ocean  surface  on  a  daily  basis.  Thus, 
surface  safellife  observations  and  an  effecfive  mefhod 
fo  assimilafe  fhese  dafa  are  a  requiremenf. 

In  fhe  pasf  few  years,  compufer  fechnology  has 
developed  fo  fhe  poinf  where  a  sysfem  capable  of 
resolving  fhe  oceanic  mesoscale  (i.e.  eddy-resolving) 
can  be  run  on  a  daily  basis.  This  firsf  generafion  global 
ocean  prediction  capabilify  has  been  developed  af  fhe 
Naval  Research  Laborafory  (NRL)  and  all  componenfs 
of  fhe  sysfem  are  rurming  in  real-time  af  fhe  Naval 


Oceanographic  Office  (NAVOCEANO).  NRL  research 
dating  back  fo  fhe  early  1980s,  which  focused  on  high 
horizonfal  resolution,  safellife  alfimefer  dafa  assimila- 
fion,  and  sfafisfical  downward  projection  of  surface 
safellife  observafions  has  led  fo  a  sysfem  made  up  of 
fhree  main  componenfs:  fhe  Modular  Ocean  Dafa 
Assimilafion  Sysfem  (MODAS),  fhe  NRL  Layered 
Ocean  Model  (NLOM),  and  fhe  NRL  Coasfal  Ocean 
Model  (NCOM).  This  sysfem  provides  fhe  Navy  wifh 
fhe  firsf  global  ocean  analysis /prediction  capabilify  for 
fleef  supporf  and  is  also  a  confribufion  fo  fhe  opera- 
fional  phase  of  fhe  multi-national  Global  Ocean  Data 
Assimilation  Experiment  (GODAE;  International 
GODAE  Steering  Team,  2000). 

The  delivery  of  fhe  firsf  generafion  sysfem  began  in 
1999  wifh  a  global  ocean  analysis  capabilify  (MODAS), 
followed  in  2000  wifh  fhe  fransifion  of  a  Yu  horizonfal 
resolution  global  ocean  model  (NLOM)  fhaf  became  fhe 
firsf  operational  sysfem  fo  properly  resolve  mesoscale 
ocean  d5mamics.  If  will  be  complefed  by  fhe  delivery  of 
a  /4°  global  version  of  fhe  NCOM,  wifh  higher  verfical 
resolution  in  fhe  mixed  layer,  for  improved  upper- 
ocean  predicfion  and  boundary  condifions  for  higher 
resolution  coasfal  models  (Eigure  1). 

There  is  sfrong  evidence  fhaf  ocean  models  need  fo 
use  grid  cells  for  each  predicfed  variable  fhaf  are  af 
mosf  abouf  7  km  across  af  mid-lafifudes.  NRL  research 
has  shown  fhaf  doubling  fhe  horizonfal  resolution  fo 
3.5  km  per  cell  gives  subsfanfial  improvemenf  buf  dou¬ 
bling  again  fo  1.7  km  gives  only  modes!  additional 
improvemenf  (Hurlburf  and  Hogan,  2000).  Eor  fhe 
NLOM  grid,  fhese  resolufions  franslafe  fo  1^6°,  Y^-i  and 
1^4°,  respectively.  This  is  for  fhe  global  and  basin-scale. 
Much  finer  resolufion  is  required  for  coasfal  regions.  Af 
3.5  km,  fhe  optimal  resolution  is  finer  fhan  mighf  be 


Oceanography  •  Vol.  15  •  No.  1/2002 


29 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

2QQ  j  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2001  to  00-00-2001 

4.  TITLE  AND  SUBTITLE 

Navy  Real-time  Global  Modeling  Systems 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Oceanographic  Division, Stennis  Space 

Center, MS, 39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

15 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Figure  1.  Snapshots 
of  sea  surface  height  from 
Vw’  real-time  NLOM 
(upper  panel)  and  sea 
surface  temperature  from 
V/  NCOM  (lower 
panel)  for  May  19,  2001. 
Note  that  NLOM’s  pole- 
ward  extent  is  72°S  and 
65°N,  while  NCOM  is 
fully  global  including  the 
Arctic. 


expected  based  on  the  size  of  eddies.  In  relation  to 
ocean  eddy  size  it  is  similar  to  the  resolution  currently 
used  by  the  leading  weather  forecasting  models  in  rela¬ 
tion  to  the  size  of  afmospheric  highs  and  lows. 

More  specifically,  our  research  has  shown  fhaf  fine 
resolution  of  fhe  ocean  eddy  scale  is  required  fo  obfain 
dynamic  coupling  between  upper  ocean  currenfs  and 
seafloor  fopography  via  flow  insfabilifies,  a  mecha¬ 
nism  fhaf  occurs  wifhouf  direcf  confacf  between  fhe 
upper  ocean  currenfs  and  fhe  fopography.  This  cou¬ 
pling  can  sfrongly  affecf  fhe  pafhways  of  upper  ocean 
currenfs  and  fronfs,  including  fhe  Gulf  Sfream  in  fhe 
Aflanfic,  fhe  Kuroshio  in  fhe  Pacific  (Hurlburf  ef  al., 
1996;  Hurlburf  and  Mefzger,  1998),  fhe  Easf  Korea 
Warm  Currenf  in  fhe  Japan/Easf  Sea  (Hogan  and 
Hurlburf,  2000)  and  fhe  Tasman  Eronf  befween 
Ausfralia  and  New  Zealand  (Tilburg  ef  al.,  2001).  The 
high-resolufion  is  also  required  fo  obfain  sharp  fronfs 
fhaf  span  major  ocean  basins  (Hurlburf  ef  al.,  1996)  and 
for  adequafe  represenfafion  of  sfraifs  and  islands  (e.g. 
Mefzger  and  Hurlburf,  2001).  In  addition,  if  can  affecf 
fhe  large-scale  shape  of  ocean  gyres  such  as  fhe 
Sargasso  Sea  in  fhe  Aflanfic  (Hurlburf  and  Hogan, 
2000).  As  discussed  lafer,  simulation  skill  for  ocean 
mesoscale  feafures  is  crucial  for  an  ocean  model  fo  acf 
as  an  effective  dynamic  inferpolafor  for  safellife  alfime- 


fer  dafa  and  for  successful  forecasfing  of  ocean  eddies 
and  fhe  meandering  of  ocean  currenfs  and  fronfs. 

NRL  Global  Ocean  Modeling  System 

Motivation  for  NLOM-NCOM  System 

NRL's  global  ocean  prediction  strategy,  as  dis¬ 
cussed  by  Harding  et  al.  (1999),  is  driven  by  the  need 
for  ocean  models  wifh  fhe  high  horizonfal  resolution 
required  fo  successfully  simulafe  fhe  variabilify  of 
mesoscale  feafures  and  fhe  high  vertical  resolution 
required  near  fhe  surface  fo  resolve  fhe  physics  of  fhe 
upper  ocean.  On  global  scales  if  is  nof  yef  possible  fo 
run  affordably  one  ocean  model  in  real  time  fhaf  has 
fhe  horizonfal  and  vertical  resolufion  needed  fo  meef 
fhese  requiremenfs.  Thus,  fhe  firsf  generafion  global 
ocean  prediction  sysfem  uses  a  fwo-model  approach 
based  on  NLOM  and  NCOM. 

NLOM,  forced  by  daily  oufpuf  from  fhe  Navy 
Operational  Global  Afmospheric  Predicfion  Sysfem 
(NOGAPS;  see  Rosmond  ef  al.,  fhis  issue),  assimilafes 
alfimefer-derived  sea  surface  heighf  (SSH)  wifh  fhe 
high  horizonfal  resolufion  needed  for  fhe  mesoscale 
d5mamics  and  serves  as  an  effecfive  dynamic  inferpola¬ 
for  of  fhe  alfimefer  dafa  in  mapping  fhis  variabilify. 
NCOM,  also  forced  by  NOGAPS  and  deriving  ifs 
mesoscale  information  from  fhe  3-dimensional 
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Figure  2.  Wiring  diagram  ofNRL  global  ocean  prediction  system  baseline  showing  relationship  between  MOD  AS,  NLOM, 
and  NCOM  modules  and  the  role  of  each  component  in  the  overall  system. 


MODAS  temperature  and  salinity  profiles  derived 
from  fhe  NLOM  SSH  and  fhe  MODAS  SST  analyses, 
provides  a  global  oceanic  boundary  layer  predicfion 
capabilify  from  fhe  open  ocean  and  info  fhe  coasf.  Wifh 
a  sophisficafed  surface  mixed-layer  formulafion 
(Mellor  and  Yamada,  1982)  and  fhe  abilify  fo  run  a  fully 
global  grid  including  fhe  Arcfic  and  confinenfal 
shelves,  NCOM  can  be  used  fo  provide  boundary  and 
inifial  conditions  fo  nesfed  higher  resolution  coasfal 
models  and  can  also  be  used  for  global  air/  sea  inferac- 
fion  applications.  Until  NCOM,  or  ofher  similar  mod¬ 
els,  can  run  on  operational  compufers  af  sufficienf  hor- 
izonfal  resolution  fo  resolve  fhe  ocean  mesoscale,  we 
expecf  fhe  fwo-model  approach  will  be  necessary. 

Figure  2  shows  a  schematic  of  fhe  enfire  global  sys- 
fem  emphasizing  fhe  individual  roles  of  NLOM, 
MODAS,  and  NCOM.  The  analysis  portion  of  fhe  sys- 
fem  is  based  on  MODAS  (Fox  ef  al.,  2001),  which  is 
now  operational  af  NAVOCEANO.  MODAS  is  an  opti¬ 
mum  inferpolafion  analysis  sysfem  wifh  fhe  abilify  fo 


generafe  S5mfhefic  femperafure  and  salinify  profiles 
based  on  observed  sea  surface  femperafure  (SST)  and 
SSH.  MODAS  confains  a  global  bimonfhly  dafabase  of 
regression  coefficienfs  relating  deviations  of  SSH  and 
SST  from  climafology  fo  deviafions  of  subsurface  fem¬ 
perafure  from  climafological  profiles.  This  S5mfhefic 
profile  capabilify  is  parf  of  a  global  climafology  fhaf 
exfends  info  shallow  wafer  using  only  SST  where  fhe 
depfh  is  <  200  mefers.  The  currenf  operational  global 
MODAS  capabilify  generafes  synfhefic  profiles  based 
on  global  SSH  and  SST  analyses  performed  daily  af 
fhe  NAVOCEANO  Alfimefry  Dafa  Eusion  Cenfer 
(ADEC;  see  Jacobs  ef  al.,  fhis  issue).  The  global  MODAS 
capabilify  will  be  upgraded  fo  perform  a  full  3-D  analy¬ 
sis  by  assimilating  available  in  situ  dafa  and  fhe  analy¬ 
sis  profiles  will  be  used  for  assimilafion  info  fhe  global 
NCOM  (discussed  lafer). 

Operational  applications  for  fhe  NLOM  sysfem 
include  improved  nowcasfs  (i.e.  fhe  currenf  sfafe  of  fhe 
ocean)  and  forecasfs  of  mesoscale  variabilify  including 
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major  front  and  eddy  systems  like  the  Gulf  Stream  and 
Kuroshio  Extension.  These  will  be  used  by  NAVO- 
CEANO  to  improve  their  capability  to  depict  the  loca¬ 
tions  of  oceanic  fronts  and  eddies.  Also,  accurate  high- 
resolution  NLOM  SSH  products  will  improve  MOD  AS 
3-D  temperature  and  salinity  synthetic  profiles,  thus 
providing  a  more  accurate  and  dynamically  realistic 
initialization  for  global  NCOM.  NLOM  provides  a 
"first  of  its  kind"  forecast  capability  for  ocean  fronts 
and  eddies,  including  the  ability  to  generate  up  to  30- 
day  forecast  MOD  AS  S5mthetic  profiles  generated  from 
NLOM's  high-resolution  SSH  and  SST  fields. 

NCOM  will  provide  an  ocean  nowcast /forecast 
capability  for  up  to  5-days  that  will  be  used  as  bound¬ 
ary  conditions  for  nested  high-resolution  fixed  and 
relocatable  coastal  forecasting  systems.  Also,  NCOM's 
more  sophisticated  surface  mixed  layer  formulation 
should  allow  the  cycling  of  NCOM  forecasts  of  SST 
and  mixed-layer  depth  with  the  MOD  AS  SST  analysis 
to  improve  the  upper  ocean  nowcasts.  Surface  current 
nowcasts  and  forecasts  from  NCOM  will  be  used  for 
search  and  rescue  operations  and  also  for  optimum 
track  ship  routing.  Studies  are  currently  imderway  at 
NRL  evaluating  NCOM  performance  for  upper  ocean 
nowcast  and  forecast  skill.  The  skill  assessment  will 
provide  measures  of  performance  and  model  error  that 
can  be  used  to  improve  the  real-time  assimilation  and 
provide  a  baseline  for  future  upgrades  to  the  global 
modeling  system. 

NRL  Layered  Ocean  Model  System 

The  NRL  Layered  Ocean  Model  is  a  descendent  of 
the  primitive  equation  model  by  Hurlburt  and 
Thompson  (1980).  Since  then  the  capability  of  the 
model  has  been  greatly  expanded  (Wallcraft,  1991; 
Wallcraft  and  Moore,  1997;  Moore  and  Wallcraft,  1998). 

The  model  has  a  free  surface,  realistic  topography 
and  a  nearly  global  domain  that  extends  from  72°S  to 
65°N  for  all  deep  water  (>  200  m)  regions  and  selected 
straits  (e.g.  Tsushima).  The  horizontal  resolution  of 
each  model  variable  is  Yu  in  latitude  by  “7512°  in  longi¬ 
tude  and  at  this  resolution,  the  model  is  eddy-resolv¬ 
ing.  The  model  is  thermod5mamic,  meaning  it  contains 
a  horizontally  and  temporally  varying  density  variable 
that  can  be  modified  by  known  physical  processes  (i.e. 
atmospheric  thermal  forcing,  advection,  diffusion, 
entrainment)  or  by  relaxation  to  a  density  climatology 
such  as  MODAS.  The  latter  helps  NLOM  maintain 
accurate  SSH.  Unlike  models  with  fixed  levels  in  the 
vertical,  NLOM's  density  relaxation  does  not  damp  the 
SSH  anomalies  because  most  of  the  information  about 
these  anomalies  is  carried  in  NLOM's  layer  thickness 
variations.  Eor  example,  NLOM  maintained  a  Rossby 
wave  generated  by  the  1982-83  El  Nino  for  at  least  a 
decade  (Jacobs  et  al.,  1994)  without  oceanic  data  assim¬ 
ilation  except  for  relaxation  to  climatological  density. 
The  model  equations  and  additional  details  concerning 
NLOM  can  be  found  in  Hurlburt  et  al.  (1996),  Metzger 
and  Hurlburt  (1996),  Moore  and  Wallcraft  (1998)  and 


Shriver  and  Hurlburt  (1997).  A  bulk  mixed  layer  model 
and  sea  surface  temperature  was  added  to  NLOM, 
Wallcraft  et  al.  (2002),  and  is  part  of  the  system  transi¬ 
tioned  to  NAVOCEANO. 

The  model  was  spun  up  over  a  10  year  period  from 
1990-2000.  During  the  spin  up  a  hybrid  wind  was  used 
consisting  of  a  combination  of  Eleet  Numerical 
Meteorology  and  Oceanography  Center  (ENMOC) 
NOGAPS  and  Hellerman  and  Rosenstein  (HR,1983) 
surface  stresses  where  the  long-term  mean  (August 
1990-July  1999  for  ENMOC)  is  subtracted  from  the 
ENMOC  stresses  and  replaced  with  the  HR  armual 
mean.  European  Centre  for  Medium-Range  Weather 
Eorecasts  thermal  forcing  was  used  from  1990-1997, 
followed  by  a  switch  to  NOGAPS  thermal  forcing  in 
1998.  In  1997  (1998),  SSH  (SST)  assimilation  began.  This 
was  continued  until  October  2000,  when  the  system 
was  transitioned  to  NAVOCEANO  where  it  has  been 
rurming  in  real  time  since  18  October  2000  and  became 
an  operational  product  on  27  September  2001.  It  is 
updated  daily  as  part  of  the  operational  rim  stream 
with  SSH  and  SST  assimilation  and  is  forced  by  real¬ 
time  NOGAPS  wind  and  thermal  forcing. 

NLOM  Data  Assimilation 

The  NLOM  system  assimilates  both  SSH  and  SST 
but  using  different  methods.  The  SSH  assimilation  is 
performed  daily  and  is  an  optimum  interpolation  (OI) 
based  data  assimilation  scheme.  The  technique  uses  an 
OI  deviation  analysis  with  the  model  as  the  first  guess, 
the  statistical  inference  technique  of  Hurlburt  et  al. 
(1990)  for  vertical  projection,  geostrophic  balance  for 
velocity  deviations  outside  an  equatorial  band,  and 
incremental  insertion  to  further  reduce  inertia-gravity 
wave  generation  (Smedstad  and  Pox,  1994).  Jacobs  et  al. 
(1999)  calculated  the  spatially  varying,  anisotropic 
mesoscale  covariance  functions  for  the  OI  analyses 
from  TOPEX/POSEIDON  and  ERS-2  data  over 
1993-99.  The  SST  assimilation  is  performed  by  relaxing 
NLOM  SST  to  the  real  time  MODAS  global  SST  analy¬ 
sis  performed  daily  at  NAVOCEANO. 

The  assimilation  cycle  for  the  operational  system 
goes  back  three  days  in  time  and  uses  analysis  wind 
and  thermal  forcing  as  SSH  and  SST  data  are  assimilat¬ 
ed  up  to  the  nowcast  time.  In  the  process  the  steps  for 
SSH  in  the  paragraph  above  are  repeated  once  a  day. 
This  allows  the  use  of  newly  received  data  and  more 
accurate  altimeter  data  that  get  updated,  as  more  pre¬ 
cise  orbit  calculations  are  available  to  the  ADPC.  This 
approach  improves  the  accuracy  of  nowcasts  used  to 
initialize  model  forecasts.  The  system  performs  a  daily 
4-day  forecast,  except  for  once  a  week  when  a  30-day 
forecast  is  made.  During  the  forecast  period,  SST  is 
relaxed  to  climatologically  corrected  persistence  of  the 
nowcast  SST  with  a  relaxation  time  scale  of  Yi  the  fore¬ 
cast  length  (i.e.  1  day  for  a  4-day  forecast  and  1  week  for 
a  4- week  forecast).  This  is  a  necessary  step  for  the  long¬ 
term  (30-day)  SST  forecasts  because  the  forecast  atmos¬ 
pheric  thermal  forcing  is  only  available  out  to  five  days. 
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Research  at  NRL  has  shown  this  method  produces  a 
more  accurate  SST  forecast.  During  a  30-day  forecast, 
the  position  of  oceanic  fronfs,  currenf/fronfal  mean¬ 
ders  and  eddies  can  change  subsfanfially.  The  SST 
relaxation  is  weak  enough  fhaf  if  allows  fhe  model  fo 
keep  fhe  SST  forecasfs  of  fhese  feafures  in  beffer  phase 
wifh  fhe  forecasfs  of  SSH  and  surface  currenfs  fhan  cli- 
mafologically  correcfed  persisfence.  There  is  no  relax¬ 
ation  for  SSH  because  fhe  model  is  used  fo  forecasf  fhe 
movemenf  of  fronfs  and  eddies  due  fo  flow  insfabilifies 
fhaf  are  more  sensifive  fo  initial  conditions  fhan  fo 
afmospheric  forcing. 

NRL  Coastal  Ocean  Model  System 

NCOM  Version  1.0  (Martin,  2000)  is  based  primari¬ 
ly  on  two  earlier  ocean  models,  the  Princeton  Ocean 
Model  (POM;  Blumberg  and  Mellor,  1983;  Blumberg 
and  Mellor,  1987)  and  the  Sigma/Z-level  Model  (SZM; 
Martin  et  al.,  1998).  SZM  is  similar  to  POM,  but  differs 
in  several  ways;  mosf  nofably  if  uses  a  hybrid  vertical 
coordinafe  sysfem.  The  vertical  grid  is  sef  up  fo  offer  a 
choice  of  sigma  layers  or  z-levels,  or  some  combination 
wifh  sigma  layers  in  fhe  shallow  wafer  and  z-levels  in 
fhe  deeper  wafer. 

Like  NLOM,  NCOM  is  designed  fo  be  a  scalable, 
porfable  code  fhaf  runs  efficienfly  on  differenf  compu- 
fafional  archifecfures.  The  model  was  designed  fo  fif 
wifhin  fhe  very  specific  code  sfrucfure  of  fhe  Coupled 
Ocean  Afmosphere  Mesoscale  Prediction  Sysfem 
(COAMPS)  fhaf  provides  for  calls  fo  (i.e.,  coupling 
between)  bofh  fhe  afmospheric  and  ocean  models 
wifhin  fhe  same  Forfran  program  (see  Hodur  ef  al.  in 
fhis  issue).  Thus,  NCOM  is  already  sef  up  for  modeling 
coupled  air-ocean  sysfems.  NCOM  also  provides  for  an 
arbifrary  number  of  nesfing  levels  wifhin  fhe  same 
Forfran  program. 

The  fransifion  of  global  NCOM  fo  NAVOCEANO 
is  planned  by  fhe  end  of  2002.  This  version  is  on  a 
curvilinear  global  grid  exfending  from  80°S  fo  a  com- 
plefe  Arcfic  cap.  This  cap  is  a  bipolar  grid  wifh  poles 
over  land  af  47°N  in  Norfh  America  and  Asia.  Global 
NCOM  uses  realistic  fopography  including  all  confi- 
nenfal  shelves  fo  depfhs  >  5  m.  The  plarmed  opera¬ 
tional  grid  resolution  for  global  NCOM  is  -Ve"  in  longi- 
fude  by  -Ve"  X  cos  0  in  lafifude  (where  0  is  lafifude) 
excepf  for  fhe  Arcfic  cap,  which  mafches  fhe  longifudes 
of  fhe  spherical  grid  af  47°N.  This  franslafes  fo  -15-16 
km  grid  spacing  af  mid-lafifudes  or  ~Vs°  in  lafifude 
near  45°N,  fhe  nominal  mid-lafifude  grid  resolufion  in 
degrees. 

Currenfly  fhere  are  two  versions  of  global  NCOM 
fhaf  have  been  fesfed  by  NRL  and  resulfs  from  bofh 
will  be  shown.  The  firsf  version  is  a  Vr"  (af  45°N)  global 
model  wifh  20  sigma-z  levels  (6  sigma,  14  z).  If  has 
been  run  in  real-time  dafa-assimilafive  mode  since 
Augusf  2000.  This  model  is  forced  using  NOGAPS  3- 
hourly  wind  sfresses  and  heaf  fluxes  and  assimilafes 
fhe  MOD  AS  global  SST  analysis  and  3-D  femperafure 
and  salinify  S5mfhefic  profiles  derived  from  fhe  global 


MODAS  SSH  and  SST  analyses.  Currenfly  fhis  version 
produces  a  2-day  forecasf  daily  as  a  fesfbed  for  fhe 
planned  5-day  operational  forecasf.  The  Vs"  (af  45°N) 
version  has  40  sigma-z  levels  (19  sigma,  21  z)  and  has 
been  rurming  in  real-fime  dafa-assimilafive  mode  since 
October  2001  forced  wifh  NOGAPS  3-hourly  wind 
sfresses  and  heaf  fluxes.  The  Vs"  model  assimilafes  SST 
from  MODAS  and  has  been  fesfed  wifh  assimilafion  of 
synfhefic  femperafure  and  salinify  derived  from  bofh 
fhe  global  MODAS  SSH  analysis  and  fhe  global 
NLOM  SSH  nowcasf. 

NCOM  Data  Assimilation 

The  data  assimilation  for  fhe  Vs"  global  NCOM  is 
accomplished  via  two  mechanisms:  adjustment  of  sur¬ 
face  heat  and  freshwater  fluxes  and  insertion  of  sub¬ 
surface  temperature  and  salinity  profiles.  In  either  case, 
the  strength  of  assimilation  is  controlled  by  a  gridded 
weighting  function  that  reflects  the  relative  confidence 
between  the  model  and  data.  The  slow  data  insertion 
technique  allows  the  model  to  incrementally  adjust  to 
the  data  with  minimal  d5mamic  disruption.  Preparation 
of  data  fields  is  independent  of  the  NCOM  assimilation 
itself,  allowing  the  model  to  accommodate  a  variety  of 
approaches  that  prepare  the  observational  analyses. 

Present  data  preparation  for  Vs"  global  NCOM 
begins  with  SST  fields  from  MODAS  2-D  (the  opera¬ 
tional  Ys"  global  MODAS  SST  analysis  performed  daily 
at  NAVOCEANO)  and  SSH  from  Yis"  global  NLOM. 
Note  that  MODAS  2-D  also  produces  a  Vs"  SSH  analysis 
from  altimeter  data,  but  in  this  global  model  applica¬ 
tion  the  NLOM  SSH  field  is  used  to  take  advantage  of 
the  improved  nowcast  skill  of  the  NLOM  dynamics 
over  the  MODAS  statistical  model.  Subtracting  the 
NLOM  climatological  mean  SSH  from  the  daily  SSH 
yields  a  SSH  deviation,  which  is  passed  along  with  the 
MODAS  2-D  SST  to  the  MODAS  module  for  calculation 
of  synthetic  profiles  of  temperature  and  salinity.  SSH 
does  not  presently  provide  value-added  to  the  MODAS 
synthetic  profiles  in  shallow  water.  Therefore,  a  mask, 
which  combines  a  smoothed  topography  with  con¬ 
straints  imposed  by  the  boundaries  of  the  NLOM 
domain,  determines  the  extent  to  which  SSH  is  used  in 
estimating  subsurface  temperature.  Combined  SSH 
and  SST  regressions  are  used  fully  where  the  mask 
depth  is  greater  than  600  m,  while  MODAS  SST-only 
regressions  are  used  in  regions  with  mask  depth  less 
than  200  m,  not  covered  by  the  NLOM  domain.  Profiles 
are  based  on  smoothly  blended  regressions  for  regions 
of  intermediate  mask  depth. 

Salinity  is  estimated  from  temperature  according  to 
local  climatological  T-S  regressions.  At  this  stage  avail¬ 
able  in  situ  observations  may  be  assimilated  into  the  3- 
D  analyses.  After  converting  temperature  to  potential 
temperature,  the  3-D  potential  temperature  and  salini¬ 
ty  fields  are  interpolated  to  the  model  grid,  and  masked 
according  to  the  model  land  and  topographic  depth. 
Surface  salinity  is  extracted  and  paired  with  the  inter¬ 
polated  SST  for  adjusting  the  surface  fluxes. 
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In  an  operational  model  run  for  nowcast  time  T,  Ys" 
global  NCOM  is  initialized  with  the  previous  analysis 
for  time  T-3  days  to  allow  the  analysis  to  take  advan¬ 
tage  of  the  more  accurate  SSH  analysis  from  NLOM, 
which  does  a  daily  re-analysis  for  time  T-3  to  T.  Surface 
heat  and  salinity  fluxes  are  adjusted  by  the  data-model 
differences  times  a  flux  rate.  Subsurface  temperature 
and  salinity  are  nudged  toward  the  data  according  to  a 
3-D  weighting  factor. 

Real-time  results 
NLOM 

NLOM  SSH  nowcast  comparisons  with  NAVO- 
CEANO  frontal  analyses 

The  Vi6°  global  NLOM  has  been  rurming  in  real 
time  with  SSH  and  SST  assimilation  since  18  October 
2000  and  became  an  operational  product  of  NAVO- 
CEANO  on  27  September  2001.  NRL  has  developed 
evaluation  software  and  has  been  monitoring  the  per¬ 
formance  of  the  system  to  establish  the  baseline  per¬ 
formance  measurements  for  the  first  generation  opera¬ 
tional  system.  At  the  present  time  only  satellite  data  are 
assimilated,  which  allows  comparison  with  independ¬ 
ent  in  situ  data  from  buoys  and  expendable  bathyther¬ 
mographs  (XBTs).  In  the  near  future,  in  situ  data  will  be 
assimilated  into  the  MODAS/NLOM/NCOM  system, 
but  for  now  will  be  used  to  determine  the  baseline  per¬ 
formance  of  the  entire  system.  In  real-time,  the  wealth 
of  the  ocean  data  available  comes  from  global  satellites 
and  these  data  will  be  the  major  contributors  to  model 
analysis  and  forecast  skill,  so  analysis  of  the  system 
performance  based  only  on  satellite  data  is  an  impor¬ 
tant,  but  also  very  demanding  evaluation. 

One  evaluation  of  the  NLOM  system  is  its  ability 
to  nowcast  the  positions  of  major  fronts  and  eddies  on 
the  global  scale.  The  Warfighting  Support  Center 
(WSC)  at  NAVOCEANO  relies  on  satellite  Infra-Red 
(IR)  SST  data  to  locate  fronts  and  eddies  for  the  global 
ocean  and  releases  frontal  analysis  products  to  the 
fleet.  The  NLOM  system  will  allow  the  WSC  analysis 
to  use  daily  nowcasts  and  animations  of  SSH  as  anoth¬ 
er  powerful  guidance  tool  to  improve  the  quality  of 
frontal  analysis  products.  This  is  particularly  signifi¬ 
cant  because  SSH  is  a  better  indicator  of  subsurface 
frontal  location  than  SST.  An  example  where  this  has 
an  impact  is  discussed  later  in  this  subsection. 
Specifically,  NLOM  provides  a  daily  map  of  the  ocean 
mesoscale  SSH  field  that  can  be  used  as  a  tool  to  help 
the  WSC  interpret  cloud  filled  IR  images.  In  addition, 
by  using  animations  of  the  NLOM  SSH  field,  the  ana¬ 
lyst  can  better  track  front  and  eddy  movements  to  help 
analyze  the  space  and  time  continuity  of  the  ocean 
mesoscale  in  areas  where  frontal  analysis  is  required. 

Eigure  3  shows  the  NLOM  SSH  nowcast  for  the 
Gulf  Stream  overlaid  with  a  portion  of  the  daily  WSC 
frontal  analysis  product  for  Jime  4, 2001  (Eigure  3a)  and 
for  June  11,  2001  (Eigure  3b).  One  must  be  careful  in 
interpreting  these  results  because  of  the  imcertain 
quality  of  the  WSC  frontal  analysis.  Due  to  clouds,  the 


analyst  may  not  have  seen  the  front  for  several  days 
and  may  be  depending  on  a  previous  (older)  frontal 
analysis  or  may  be  analyzing  for  a  shallow  surface  fea¬ 
ture  of  SST  off  the  Gulf  Stream.  In  this  case  for  June,  the 
WSG  had  good  quality  images  and  the  frontal  position 
has  high  confidence  along  most  of  the  Gulf  Stream 
path.  However,  the  WSG  analysis  did  not  use  the 
NLOM  model  output  and  can  be  considered  a  com¬ 
pletely  independent  analysis. 

Eor  both  June  4  and  June  11,  2001,  the  NLOM  Gulf 
Stream  pathway  has  good  agreement  with  the  WSG 
frontal  analysis.  Both  NLOM  and  the  frontal  analysis 
agree  on  the  position  and  shape  of  the  meanders  at 
67°W  and  53°W,  features  that  change  their  shape 
between  June  4  and  June  11.  During  this  period  the 
meanders  deepen,  as  seen  in  both  the  model  and  the 
frontal  analyses.  The  model  also  captures  a  warm  ring 
shedding  event  near  62°W  that  is  verified  by  the  WSG 
frontal  analysis.  On  June  4,  both  products  indicate  a 
deep  ridge  near  60°W  that  bends  back  toward  the  west 
and  on  June  11,  both  indicate  that  a  warm  eddy  has 
shed  from  the  Gulf  Stream.  There  is  also  good  agree¬ 
ment  on  a  trough  at  47°W,  which  indicates  that  the 
model  has  a  realistic  pathway  from  Gape  Hatteras  to 
the  Grand  Banks  when  compared  with  the  WSG  frontal 
analysis.  Additional  comparisons  are  available  on  the 
NLOM  web  site  at  http://www7320.nrlssc.navy.mil- 
/  global_nlom. 

The  real  value  of  the  NLOM  SSH  frontal  analysis 
will  be  during  times  when  the  IR  data  is  cloudy  or 
when  a  subsurface  front  is  masked  by  uniform  SST  as 
occurs  for  the  Gulf  of  Mexico  Loop  Gurrent  in  summer, 
for  example.  The  WSG  plans  to  use  the  NLOM  product 
as  a  first  guess  for  their  frontal  analysis  product.  The 
errors  in  the  WSG  analyses  are  mainly  due  to  IR 
imagery  that  is  difficult  or  impossible  to  interpret 
because  of  clouds  or  to  lack  of  a  frontal  signature  in 
SST.  NLOM  provides  a  useful  guidance  product  for  the 
WSG  analysts  to  improve  the  frontal  analysis  in  areas 
where  it  is  cloudy  or  during  the  warmer  months  when 
the  fronts  lose  their  surface  SST  signature. 

NLOM  Sea  Surface  Height  for  deriving 
MOD  AS  synthetics 

As  discussed  earlier,  one  of  the  operational  uses  for 
the  Vi6°  global  NLOM  is  to  provide  accurate  high-reso¬ 
lution  SSH  nowcasts  with  realistic  mesoscale  features 
to  derive  MODAS  S5mthetic  temperature  and  salinity 
profiles.  The  increased  accuracy  of  the  profiles  will  help 
to  improve  the  full  MODAS  3-D  analysis,  which  will 
then  be  assimilated  into  the  global  NGOM  system  pro¬ 
viding  an  improved  initialization  for  the  NGOM  fore¬ 
casts.  Since  NLOM  has  the  needed  resolution  to  act  as 
a  dynamic  interpolator  of  the  satellite  altimeter  data,  it 
provides  the  global  NGOM  with  a  more  d5mamically 
correct  initialization  than  MODAS.  Some  of  the  advan¬ 
tages  of  the  higher  resolution  NLOM  S5mthetics  as 
input  to  global  NGOM  include,  better  defined  ocean 
fronts,  better  definition  of  boundary  currents,  western 
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Figure  3.  SSH  analysis  (nowcast)  in  the  Gulf  Stream  region  from  the  real-time  Yu  global  NLOMfor  a)  June  4,  2001  and  b) 
June  11,  2001.  Superimposed  on  each  is  an  independent  Gulf  Stream  north  wall  frontal  analysis  determined  from  satellite  IR 
imagery  (white  lines)  by  the  Naval  Oceanographic  Office  for  the  same  days.  The  color  palette  was  chosen  to  emphasize  the  loca¬ 
tion  of  the  Gulf  Stream  and  associated  eddies. 


boundary  currents  that  more  realistically  penetrate 
into  the  basin,  more  realistic  SSH  gradients  across 
fronts,  and  more  realistic  shapes  and  size  of  eddies. 

Figure  4  shows  the  results  of  comparisons  between 
MOD  AS  and  NLOM  generated  3-D  S5mthetics  vs.  inde¬ 
pendent  XBTs  from  the  NLOM  Operational  Test  per¬ 
formed  at  NAVOCEANO.  A  total  of  about  4900  XBTs 
and  conductivity-temperature-depth  (CTD)  measure¬ 
ments  collected  worldwide  in  January  and  February 
2001  were  used  for  independent  evaluation  of  the  syn¬ 
thetic  profiles  calculated  using  the  MODAS  2-D  SSH 
field  and  the  NLOM  SSH  field.  None  of  the  measured 
profiles  were  assimilated  into  either  MODAS  2-D  or 
NLOM.  The  test  was  to  compare  the  root  mean  square 


(RMS)  error  (°C)  of  S5mthetic  temperafure  profiles  cal- 
culafed  from  MODAS  and  NLOM  SSH  (vs.  XBT  fem- 
perature  profiles).  Errors  were  calculated  as  the  differ¬ 
ence  between  a  measured  profile  and  ifs  nearest 
synthetic  neighbor  at  a  series  of  sfandard  depfhs.  This 
fest  was  performed  on  the  entire  data  set  (Figure  4a), 
and  again  on  a  subset  of  the  data  (Figure  4b)  that  con¬ 
sists  of  those  XBTs  (-800)  where  the  corresponding 
nowcast  SSH  (both  NLOM  and  MODAS)  was  greater 
than  7  cm  different  from  climatology.  The  subset  crite¬ 
ria  were  chosen  to  highlight  the  larger  improvement  of 
NLOM  and  MODAS  over  climatology  in  areas  where 
there  was  observed  mesoscale  variability. 

Figure  4a  shows  the  RMS  error  for  all  fhe  XBTs  vs. 
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operational  MOD  AS  2-D,  NLOM  and  climatology  (the 
latter  being  the  accepted  first  order  comparison  for 
ocean  analysis  sysfems).  The  resulfs  indicafe  fhaf 
NLOM  has  equal  or  lower  RMS  error  fhan  MODAS 
down  to  -400  meters  and  is  basically  the  same  at  the 
lower  depths,  but  both  NLOM  and  MODAS  2-D  are 


Figure  4.  RMS  temperature  error  (°C)  vs.  depth  (meters) 
for  MODAS  climatology  (black)  and  MODAS  synthet¬ 
ic  profiles  calculated  using  MODAS  2-D  SSH  (blue) 
and  NLOM  SSH  (red)  compared  with  a)  4900  unassim¬ 
ilated  global  XBTsfrom  January  and  February,  2001  and 
b)  a  subset  of  the  global  XBTs  (~800)  where  MODAS 
and  NLOM  SSH  at  the  XBT  location  were  both  >7  cen¬ 
timeters  away  from  the  climatological  SSH. 


superior  to  climatology.  Figure  4B  is  the  result  for  fhe 
subsef  of  XBTs  showing  fhe  comparison  with  SSH  for 
NLOM  and  MODAS  for  cases  where  there  was 
observed  mesoscale  variability  >7  cm  from  climafology. 
For  this  case,  the  value  added  of  NLOM  and  MODAS 
2-D  is  much  more  apparent.  Both  MODAS  and  NLOM 
synthetics  have  lower  RMS  error  than  climatology  and 
similar  to  the  previous  case  with  all  the  data,  NLOM 
has  equal  or  lower  errors  than  MODAS  2-D  down  to 
450  meters.  In  this  case,  MODAS  2-D  is  slightly  better 
from  500  to  800  meters  with  the  two  products  essential¬ 
ly  the  same  down  to  1000  meters. 

NLOM  SSH  forecasts 

The  ability  of  NLOM  to  forecast  SSH  and  the  posi¬ 
tions  of  major  fronfs  and  eddies  represenfs  a  new  Naval 
product  available  for  future  operational  plarming  and 
to  help  users  gauge  the  quality  of  fhe  producf.  The 
fufure  positions  of  major  ocean  fronfs  can  give  fhe 
warfighfer  guidance  on  how  changes  in  fhe  ocean  envi- 
ronmenf  mighf  affecf  fufure  missions.  An  accurafe  fore- 
casf  of  SSH  would  allow  the  Navy  to  predict  changes  in 
locations  of  mesoscale  features  (fronts  and  eddies)  that 
affect  the  3-D  temperature  and  salinity  field  by  using 
the  predicted  NLOM  SSH  and  SST  to  derive  MODAS 
synthetic  profiles.  Changes  in  the  3-D  ocean  environ¬ 
ment  would  change  acoustic  conditions  in  the  deep 
water  and  could  be  used  as  a  plarming  tool  for  future 
deployment  of  operafional  fleef  assefs  or  when  a  fleef 
asset  was  unable  to  receive  timely  products. 

Figure  5  shows  a  SSH  example  of  a  nowcasf  for 
January  3,  2001  (Figure  5a),  the  15-day  forecast  (Figure 
5b)  and  verifying  analysis  for  January  18,  2001  (Figure 
5c),  and  fhe  30-day  forecasf  (Figure  5d)  and  verifying 
analysis  (Figure  5e)  for  February  2,  2001  for  fhe 
Kuroshio  Extension  region.  The  NLOM  system  per¬ 
forms  a  30-day  forecasf  once  a  week  and  the  verification 
statistics  are  calculated  by  comparing  the  model  daily 
prediction  to  the  corresponding  validating  analysis  for 
that  same  day.  The  Kuroshio  is  in  "meander  mode"  near 
140°E  and  in  this  case  the  model  predicts  a  widening  of 
fhe  meander  after  30  days,  which  is  also  seen  in  the  ver¬ 
ifying  analysis.  The  model  also  correcfly  predicts  the 
growth  of  fhe  ridge  af  149°E  and  fhe  deepening  of  fhe 
trough  meander  near  147°E,  as  well  as  the  movement  of 
the  warm  and  cold  eddies  near  150°E  that  will  be  evi¬ 
dent  in  the  cross-section  (Eigure  7)  discussion. 

Eigure  6  shows  mean  30-day  forecast  verification 
statistics  for  SSH  RMS  error  and  SSH  anomaly  correla¬ 
tion  versus  forecasf  length  for  fhe  19  monthly  forecasts 
made  each  week  from  December  20,  2000  to  May  16, 
2001.  Anomaly  correlation  is  a  standard  statistical  cor¬ 
relation  of  deviafions  from  the  mean  SSH  that  is  similar 
to  that  used  in  atmospheric  forecast  verification  where 
values  >  0.6  are  considered  skillful  (Murphy  and 
Epstein,  1989).  Results  are  shown  for  fhe  entire  global 
domain  (Eigures  6a  and  6b),  fhe  Gulf  Stream  region 
(Eigures  6c  and  6d),  and  the  Kuroshio  region  (Eigures 
6e  and  6f).  These  statistics  are  calculated  by  comparing 
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Figure  5.  SSH  (cm)  for  the  Kuroshio  region  from  the  Vu 
global  NLOM  running  in  forecast  mode  for  a  30-day 
forecast.  The  model  forecast  was  initialized  on  January  3, 
2001  a)  and  run  for  30  days  with  forecast  b)  and  vali¬ 
dating  analysis  c)  for  January  18,  2001  and  forecast  d) 
and  validating  analysis  e)  for  February  2,  2001.  The 
model  forecast  was  forced  with  FNMOC  wind  and  ther¬ 
mal  forcing  out  through  five  days  and  then  gradually 
switched  to  climatological  forcing  for  the  remainder  of 
the  30-day  forecast.  Solid  black  line  superimposed  along 
150°E  on  each  is  the  location  of  3-D  temperature  cross 
section  shoivn  in  Figure  7. 


the  model  forecast  SSH  field  vs.  fhe  analysis  SSH  from 
NLOM  valid  on  fhe  day  of  fhe  forecasf.  For  bofh  RMS 
error  and  anomaly  correlafion  on  fhe  global  domain, 
fhe  model  is  healing  persisfence  and  climafology  (for 
RMS  error)  wifh  fhe  anomaly  correlafion  showing 
model  skill  ouf  beyond  30  days. 

Verificafion  sfafisfics  for  fhe  Gulf  Sfream  and 
Kuroshio  areas  indicafe  fhe  skill  of  fhe  model  in  regions 
dominafed  by  flow  insfabilifies  and  high  SSH  variabil- 
ify  In  fhe  Gulf  Sfream  area,  NLOM  shows  forecasf  skill 
fhaf  is  less  fhan  fhe  global  case.  In  ferms  of  RMS  SSH 
error,  NLOM  heals  persisfence  in  fhe  Gulf  Sfream  ouf  fo 
30  days.  However,  fhe  SSH  anomaly  correlafions  indi¬ 
cafe  fhaf  NLOM  fypically  has  skill  for  less  fhan  20  days. 
The  Kuroshio  resulfs  are  better  fhan  fhe  Gulf  Sfream 
and  more  similar  fo  fhe  global  case  wifh  skill  ouf 
beyond  30  days.  As  expecfed,  fhe  RMS  error  is  higher  in 
fhe  Kuroshio  fhan  fhe  global  domain  due  fo  fhe  high 
SSH  variabilify  fhaf  is  observed  in  fhaf  area  due  fo  flow 
insfabilifies.  Alfhough  nof  shown  here,  NLOM  SST  pre- 
dicfion  measures  are  also  calculafed  for  30-day  fore- 
casfs  and  show  fhaf  NLOM  on  average  heals  bofh  per¬ 
sisfence  and  climafology  for  fhe  global  domain  because 
of  fhe  relaxation  fo  climafologically  correcfed  persisf¬ 
ence  as  discussed  earlier.  Real-time  SST  prediction 
measures  are  similar  fo  fhose  shown  for  SSH.  In  addi¬ 
tion,  nowcasf  SST  comparisons  wifh  unassimilafed 
buoy  dafa  are  updafed  daily  and  available  on  NRL's 
public  NLOM  websife  discussed  earlier. 

The  difference  in  model  skill  befween  fhe  Kuroshio 
and  Gulf  Sfream  is  due  fo  fhe  resolution  of  fhe  NLOM 
sysfem  and  fhe  d5mamics  of  fhe  two  sfrong  wesfern 
boundary  currenfs.  As  sfafed  earlier,  fhe  resulfs  shown 
here  are  consisfenf  wifh  ofher  NRL  research  (Hurlburf 
and  Hogan,  2000),  which  indicafes  fhaf  af  leasf  ViT  is 
required  for  proper  simulation  of  fhe  Gulf  Sfream.  The 
forecasf  sfafisfics  indicafe  fhaf  fhe  Kuroshio  dynamics 
are  well  represenfed  by  fhe  model.  For  mosf  of  fhe 
ocean,  Vit  resolution  is  enough  fo  allow  skillful  fore- 
casfs  of  SSH.  Even  fhough  fhe  currenf  sysfem  does  nof 
show  30-day  skill  for  fhe  Gulf  Sfream,  fhe  resulfs  indi¬ 
cafe  fhaf  fhe  sysfem  has  skill  fo  af  leasf  18  days.  NRL  is 
currenfly  working  on  a  Vsz"  global  sysfem  ,  which  is  a 
possible  fufure  upgrade  fo  fhe  NLOM  sysfem,  and  fhis 
sysfem  should  improve  fhe  Gulf  Sfream  resulfs. 

Figure  7  shows  an  example  of  MOD  AS  S5mfhefic 
femperafure  and  salinify  calculafed  from  fhe  NLOM 
SSH  and  SST  analyses  and  forecasfs  on  January  3 
(Figure  7a),  January  18  (Figures  7b  and  7c),  and 
February  2,  2001  (Figures  7d  and  7e)  for  fhe  Kuroshio 
cross  secfion  along  150°E,  depicfed  by  fhe  solid  black 
line  in  Eigure  5.  Similar  fo  fhe  SSH  forecasf  comparison 
in  Eigure  5,  fhe  15  and  30-day  forecasfs  are  compared  fo 
fhe  corresponding  analysis  valid  on  fhe  forecasf  day. 
There  are  several  inferesfing  evenfs  fhaf  are  correcfly 
capfured  by  fhe  model  forecasf.  The  Kuroshio  fronf  is 
displaced  fo  fhe  norfh  during  fhe  forecasf  as  a  ridge 
meander  builds  along  150°E.  On  January  3,  2001, 
(Eigure  7a)  fhe  Kuroshio  fronf  is  af  35°N  and  moves 
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Figure  6.  Mean  SSH  forecast  verification  statistics  for  19  weekly  30-day  forecasts  from  December  20,  2000  to  May  16,  2001 
for  the  Yu  global  NLOM.  Left  column  shows  mean  SSH  RMS  error  (cm)  and  right  column  shows  mean  anomaly  correlation 
vs.  forecast  length  (days)  for  NLOM  forecast  (red  curve),  persistence  forecast  (blue  curve)  and  climatology  forecast  (black 
curve).  The  top  row  a,  h)  is  for  the  global  domain,  the  middle  row  c,  d)  is  the  Gulf  Stream  region  and  the  bottom  row  e,f)  is 
for  the  Kuroshio  region. 
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Figure  7.  Temperature  (°C)  vs.  depth  along  150°Efor  the  Kuroshio  region  derived  from  Yu  global  NLOM  SSH  and  SST  fields 
using  the  MOD  AS  synthetic  climatology  (statistics  derived  from  the  historical  hydrographic  database).  The  model  forecast  was 
initialized  on  January  3,  2001  a)  and  run  for  30  days  with  forecast  b)  and  validating  analysis  c)for  January  18,  2001  and  fore¬ 
cast  d)  and  validating  analysis  e)  for  February  2,  2001. 
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north  of  36°N  on  February  2,  2001,  which  is  seen  in 
both  the  analysis  (Figure  7e)  and  model  forecast 
(Figure  7d).  Also,  a  warm  ring  seen  in  the  analysis  near 
38°N  on  January  3  is  not  visible  in  the  model  forecast 
or  verifying  analyses  on  January  18  or  February  2  as 
the  eddy  has  propagated  to  the  west  during  the  fore¬ 
cast  period.  Similarly,  a  cold  eddy  that  was  not  in  the 
cross  section  on  January  3  is  evident  near  32°N  in  both 
the  forecast  and  analysis  on  January  18  and  February  2. 
This  case  study  indicates  how  the  3-D  temperature, 
and  consequently  the  derived  sound  speed  can  change 
and  it  highlights  the  value  of  a  skillful  long-range  fore¬ 
cast  able  to  predict  possible  oceanographic  and 
acoustic  changes  in  a  future  operational  area. 

NCOM 

Global  NCOM  SST  prediction  measures 

One  of  the  important  operational  uses  for  global 
NCOM  is  for  short-term  predictions  of  SST  and  mixed- 
layer  depth  with  high  vertical  resolution.  Using  a  more 


sophisticated  upper  ocean  mixed-layer  formulation, 
NCOM  is  designed  for  shorter  forecasts  (5-7  days)  and 
uses  thermal  forcing  from  the  operational  NOGAPS 
model.  This  contrasts  with  the  NLOM  SST  prediction 
capability  designed  for  long-term  (30  day)  forecasts 
and  depending  on  a  weak  relaxation  to  a  climatologi- 
cally  corrected  persistence  because  there  is  no  skillful 
long-term  atmospheric  heat  flux  prediction  capability 
out  to  30  days.  The  expected  NCOM  skill  for  short-term 
mixed  layer  prediction,  including  SST  forecasts,  also 
provides  the  potential  for  its  use  as  part  of  a  global  cou¬ 
pled  air /ocean  system  and  as  a  first  guess  to  the 
MOD  AS  analysis  in  a  cycling  analysis  /  forecast  system. 

As  a  testbed  for  the  development  of  the  Vs‘  global 
NCOM  capability,  NRL  has  run  a  real-time  Vi  global 
version  of  NCOM  since  August  2000  that  produces  a 
daily  nowcast  and  2-day  forecast.  SST  prediction  meas¬ 
ures  are  similar  to  those  used  for  NLOM  SSFl  where  the 
SST  forecast  is  compared  to  the  corresponding  analysis 
for  the  same  day.  Figure  8  shows  the  results  for  40  daily 
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Figure  8.  Mean  SST  forecast  verification  statistics  for  40  tzvo-day  forecasts  from  January  4,  2001  to  February  12,  2001  for  the 
Vi  global  NCOM.  Mean  SST  RMS  error  (°C)  for  NCOM  forecast  a)  and  persistence  b)  in  the  Gulf  Stream  region;  and  for 
NCOM  forecast  c)  and  persistence  d)  in  the  Agulhas  retroflection  region. 
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Figure  9.  SSH  (cm)  nowcast  for  the  Kuroshio  region  on  May  19,  2001  from  a)  Vs°  global  NCOM  and  h)  Vu  global  NLOM. 
NCOM  model  assimilates  MODAS  synthetic  temperature  and  salinity  profiles  derived  using  the  NLOM  SSH  field. 


forecast  cycles  from  January  4  to  February  12,  2001. 
Figure  8a  and  8b  are  the  mean  RMS  error  (°C)  for  the 
model  forecast  (A)  and  persistence  (B)  in  the  North 
Atlantic  region.  The  forecast  has  lower  errors  than  per¬ 
sistence  throughout  the  region  with  the  model  showing 
highest  skill  in  the  area  of  the  Gulf  Stream  front. 
Similarly  for  the  Agulhas  retroflection  area  (Figures  8c 
and  8d),  the  model  RMS  error  is  lower  than  persistence 
throughout  the  basin  with  highest  skill  in  areas  of 
frontal  instabilities. 

These  results  are  consistent  globally  and  for  differ¬ 
ent  seasons.  This  is  a  very  encouraging  result  for  SST 
prediction  skill  and  comparable  testing  will  be  done 
with  the  Ys"  global  NCOM.  In  addition,  the  evaluation 
software  that  has  been  developed  for  NLOM  can  be 


used  for  the  NCOM  system  to  evaluate  skill  levels  for 
short-term  SSH  and  3-D  temperature  prediction  skill. 

Vs‘  Real-time  Global  NCOM 

The  Ys"  global  NCOM  system  is  plarmed  for  transi¬ 
tion  to  NAVOCEANO  and  has  been  rurming  in  real¬ 
time  since  26  October  2001.  The  system  is  forced  by 
FNMOC  NOGAPS  and  assimilates  SST  data  as  well  as 
3-D  temperature  and  salinity  profiles  from  the  MODAS 
sjmthetic  climatology  using  the  Yu  global  NLOM  SSH 
and  the  Ys"  MODAS  2-D  SST  analyses  in  generating  the 
subsurface  S5mthetics. 

Figure  9  shows  a  result  from  the  Vs"  global  NCOM 
in  the  Kuroshio  region  for  May  19,  2001,  when  NCOM 
was  rurming  in  near  real-time.  Figure  9a  is  the  NCOM 
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SSH  nowcast  and  Figure  9b  is  the  corresponding  SSH 
nowcast  from  the  real-time  Vu  global  NLOM.  The 
NLOM  SSH  field  was  used  by  MOD  AS  to  calculate  the 
S5mthetic  temperature  and  salinity  profiles  that  were 
assimilated  into  the  global  NCOM.  The  results  indicate 
that  the  NCOM  SSH  field  matches  very  well  with 
NLOM,  i.e.  the  mesoscale  features  are  in  the  same  loca¬ 
tions  and  the  Kuroshio  frontal  path  is  similar  in  both 
models.  This  is  an  example  of  NLOM,  which  has  the 
resolution  required  to  properly  resolve  the  ocean 
mesoscale,  constraining  the  lower  resolution  NCOM  to 
allow  an  accurate  initialization  for  the  NCOM  short¬ 
term  forecast.  Having  a  more  accurate  initialization 
with  the  mesoscale  features  in  the  right  location  will 
allow  NCOM  to  provide  more  accurate  boundary  con¬ 
ditions  to  nested  higher  resolution  coastal  and  relocat¬ 
able  models.  It  will  also  help  to  improve  the  SST  fore¬ 
casts  by  improving  the  mesoscale  advection  field. 
Real-time  results  are  available  on  the  NCOM  web  site 
at  http:/ /www7320.nrlssc.navy.mil/global_ncom. 

Summary 

The  components  of  NRL's  NLOM/NCOM  global 
system  and  initial  validations  of  the  assimilative  mod¬ 
els  have  been  described  in  this  report.  The  system  con¬ 
sists  of  a  Yu  global  version  of  NLOM  that  has  been 
running  in  real-time  since  18  October  2000  and  a  Vs° 
global  version  of  NCOM  rurming  in  real-time  since  26 
October  2001.  The  NLOM  component  became  an 
operational  product  at  the  Naval  Oceanographic 
Office  on  27  September  2001.  It  includes  a  bulk 
mixed-layer  formulation  and  assimilation  based  on 
01  deviation  analysis  using  the  model  as  a  first  guess, 
subsurface  statistical  inference,  and  slow  data 
insertion.  NCOM  has  a  more  sophisticated  mixed- 
layer  model  and  assimilates  S5mthetic  temperature 
and  salinity  profiles  from  MODAS  based  on  the 
NLOM  SSH  field  and  the  MODAS  SST  analysis.  Both 
models  are  forced  by  FNMOC  wind  and  thermal  forc¬ 
ing  products  and  assimilate  SST  from  the  MODAS 
global  SST  analysis. 

The  NLOM  system  has  been  validated  for  nowcast 
and  forecast  skill  compared  to  independent  observa¬ 
tions.  Comparisons  of  MODAS  S5mthetic  temperature 
profiles  calculated  using  NLOM  SSH  and  the  global 
MODAS  2-D  analysis  SSH  were  performed.  The  com¬ 
parisons  with  independent  XBTs  showed  that  the 
NLOM  derived  S5mthetic  profiles  were  overall  better 
than  the  MODAS  derived  profiles.  Both  NLOM  and 
MODAS  2-D  derived  S5mthetics  were  superior  to  cli¬ 
matology  especially  where  there  are  significant  anom¬ 
alies,  including  mesoscale  variability. 

The  NLOM  system  was  tested  for  its  ability  to 
nowcast  locations  of  major  ocean  fronts  and  eddies  and 
will  be  used  to  improve  the  quality  of  the  NAVO- 
CEANO  frontal  analysis  products.  The  position  of 
fronts  from  the  NLOM  SSH  field  showed  good  agree¬ 
ment  with  independent  NAVOCEANO  frontal  analy¬ 
ses  for  the  Gulf  Stream.  The  errors  in  the  NAVO¬ 


CEANO  analysis  are  mainly  due  to  IR  imagery  that  is 
difficult  or  impossible  to  interpret  because  of  clouds. 
The  NLOM  provides  a  useful  guidance  product  for  the 
WSC  analyst  to  improve  the  frontal  analysis  in  areas 
where  interpretation  of  IR  imagery  is  difficult  due  to 
the  clouds  or  during  the  warmer  months  when  the 
fronts  lose  their  surface  SST  signal. 

The  abilify  of  NLOM  to  forecast  the  global  SSH 
field  was  also  tested.  NLOM  30-day  forecasts  were  ver¬ 
ified  against  later  analyses  and  results  shown  for  the 
global  domain  plus  the  Kuroshio  and  Gulf  Stream  sub- 
regions.  Results  indicated  that  for  the  global  and 
Kuroshio  domains,  NLOM  showed  skill  over  persist¬ 
ence  and  climatology  for  SSH  prediction  for  at  least  30 
days  and  in  the  Gulf  Stream  region  for  ~18  days.  This  is 
consistent  with  previous  NRL  research  that  indicates 
that  Vi6°  resolution  is  adequate  for  the  Kuroshio,  but  Vii" 
resolution  is  required  for  more  accurate  representation 
of  the  Gulf  Stream. 

Initial  validations  were  shown  for  global  NCOM, 
including  SST  prediction  skill  and  the  ability  of  NCOM 
to  nowcast  the  positions  of  major  fronts  and  eddies. 
Results  indicated  that  for  the  Vi  version  of  global 
NCOM,  that  has  been  running  in  real-time  since 
August  2000,  the  model  2-day  SST  forecast  had  lower 
errors  than  persistence  when  compared  to  the  analysis 
valid  on  the  forecast  day.  The  Vs‘  version  of  global 
NCOM  is  now  running  in  real  time  assimilating 
MODAS  synthetic  temperature  and  salinity  profiles 
derived  from  the  NLOM  SSH  nowcast  and  the  MODAS 
SST  analysis.  The  NCOM  SSH  nowcast  was  shown  to 
compare  well  with  the  positions  of  mesoscale  features 
seen  in  the  NLOM  SSH  nowcast,  which  has  verified 
skill  in  nowcasting  the  position  of  ocean  fronts  and 
eddies.  ^ 
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